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Abstract: Exclusive measurements of the production of ^-mesons in the pp — > pprj 
reaction have been carried out at excess energies of 16 and 37 MeV above threshold. The 
deviations from phase space are dominated by the proton-proton final state interaction and 
this influences particularly the energy distribution of the rj meson. However, evidence is 
also presented at the higher energy for the existence of an anisotropy in the angular distri- 
butions of the ?7-meson and also of the final proton-proton pair, probably to be associated 
with D-waves in this system interfering with the dominant 5-wave term. The sign of the 
r) angular anisotropy suggests that p-exchange is important for this reaction. 



The production of r] mesons in proton-proton collisions near threshold has been mea- 
sured in recent years by three different groups [H El Ej and a fairly consistent picture has 
emerged. In comparison to the analogous pion case [3J E], S'-wave production is more 
dominant near threshold, and this is generally ascribed to the presence of the N*(1535) Sn 
isobar, whose width overlaps the rj — p threshold and which has large branching fractions 
into both i] — p and tt — p [6 . Dalitz plots obtained for the pp — > pprj reaction at low 
energy show strong deviations from phase space due to the presence of the proton-proton 
final state interaction (FSI) In addition there are residual effects but the precision of 
these measurements was limited, in part, by the uncertainty in the determination of the 
proton angles. This defect has been overcome in the present experiment through the in- 
clusion of a tracking device covering the forward angles and operating in conjunction with 
the apparatus previously used. The combination allowed us to investigate the influence 
of higher partial waves in the angular distributions and to identify dependences on both 
the angles of the r] and of the proton-proton relative momentum, which become stronger 
with energy. This behaviour is the first indication of effects from higher partial waves 
in the pp — > pprj reaction. The shape of the rj angular variation is sensitive to the basic 
production mechanism and the data suggest that p-exchange provides a more important 
driving term here than 7r-exchange. 

The experiment was carried out using the PROMICE/WASA facility at the CELSIUS 
storage ring of the The Svedberg Laboratory at beam energies of 1296 and 1350 MeV, 
corresponding to centre-of-mass excess energies Q = 16 and 37 MeV respectively. Using a 
cluster gas jet target, integrated luminosities of about 50 nb _1 and 200 nb _1 were obtained 
for the two energies in a total of 30 hours of running, and these yielded about 300 and 750 
good ppr] events in the final sample. Details of the detector system are given in Refs. [HE], 
and only the main points are discussed here. 

The forward-going protons were measured in a detector system covering polar angles 
between 4° and 20° with respect to the beam direction. It consists of a tracking detector, 
followed by a three-layer scintillator hodoscope and a four-layer scintillator calorimeter. A 
second hodoscope is placed at the end of the detector system to register penetrating parti- 
cles. The tracker consists of two planes, each with four layers of thin-walled cylindrical drift 
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chambers, so-called "straw chambers", oriented in the vertical and horizontal directions. 
This arrangement allows the proton scattering angles to be reconstructed to a precision 
of better than 1° (FWHM). The r?'s are identified from their 27 decay channel, where the 
7's are detected in two CsI(Na) arrays situated on either side of the scattering chamber. 
Scintillator hodoscopes are placed in front of the Csl arrays to veto charged particles. The 
27 invariant mass resolution obtained at the ?7-meson mass is 20 MeV (RMS). 

In the off-line analysis, only those events with an identified 77 meson together with two 
energetic protons in the forward detector were retained and, in these kinematical 
fit was applied in order to extract the energy-momentum vectors for the particles. To avoid 
problems in accounting for losses through nuclear reactions in the detector material, only 
the proton directions were used, and this results in a fit with three constraints (3C). The 
number of background events with two uncorrelated 7's arising from different pions in 27r° 
production becomes negligible after applying a lower cut at 5% on the confidence level of 
the fit. In the subsequent analysis of our data, we exploit the fact that at low energies 
only a few amplitudes are allowed and that these lead to distributions in which the angular 
and momentum variables are intimately linked. It is only through the introduction of such 
a functional form into Monte Carlo simulations that we can draw any firm conclusions 
on the physics of the process. The detector response and acceptance calculations for the 
experiment were made with events generated from a full Monte Carlo simulation using 
GEANT3 [9J. The resulting geometrical acceptance for this type of event is about 0.5% 
and 2% at Q = 16 and 37 MeV respectively. The data were checked for internal consistency 
by applying different geometrical cuts. 

The sole production amplitude which survives at threshold corresponds to the transition 
3 P —> 1 So s, where we are using the standard 2S+1 Lj notation for the pp system, with the 
lower case letter denoting the angular momentum of the r?-meson with respect to the 
pp system. At slightly higher energies, amplitudes corresponding to the production of 
P- and .D-wave pp pairs introduce dependences on the pp angles with similar momentum 
threshold factors since the D-wave term can interfere with the threshold .S-wave amplitude. 
However, given that 77 production in the ^So p state is forbidden by selection rules, the first 
non-trivial 77 angular dependence is expected to come from the interference of the s- and 
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d-wave amplitudes. These considerations lead us to take the following simple form for the 
low energy pp — > pprj amplitude: 

M = A Ss <p\ (p • e-) + A Sd <f)\ (p • k)(k ■ Si) + A Ps & (q- ej) + A Ds 0} (p • • ej) . (1) 

The momenta of the initial proton and final i] in the overall cm. system are denoted by 
p and respectively, 2q is the relative momentum in the final two-proton system, ei (e/) 
the spin-one polarisation vector of the initial (final) proton-proton pair, and fa (0/) the 
corresponding spin-zero functions. The indices on the amplitudes denote the dominant 
final angular momentum state in a term. It should however be noted that the strict Ss 
partial wave amplitude is given by 

A Ss = A Ss + § k 2 A Sd + \q 2 A Ds . (2) 

Keeping only terms up to order k 2 or q 2 , i.e. a total of two units of angular momentum 
in the final intensity, the spin-averaged matrix element squared becomes 
1 



M 



\A Ss \ 2 + 2k 2 Re {A* Ss A Sd } cos 2 0„ + 2 q 2 Re {A* Ss A Ds } cos 2 6 PP + q 2 \A Ps \ 



(3) 

where 9 V and 6 PP are the angles that the r) and the pp relative momentum make with respect 
to the beam direction. 

In a model where the N*(1535) isobar is excited through one pion exchange JU]> it 
is expected that the |Ap s | 2 term should be smaller than the Re^A* Ss A Ds ^ of the Ss- 
Ds interference by a factor of the order of /i/4m p , where and m p are the the t] and 
proton masses respectively, and // = 2m p m v / (2m p + m v ) is the reduced mass of the final 
state. Since the Ps term has no characteristic angular dependence, it would be difficult to 
isolate a small effect as compared to a possible energy dependence of A Ss , and so any such 
contribution is ignored. 

The 1 Sq state of the final pp system is subject to a very strong FSI which is central 
to any analysis of low energy production. All the terms in Eq. (JHJ) are influenced by the 
FSI, with the exception of the P-wave contribution, and this will reduce even further the 
relative importance of lAp^l 2 . Without knowing the radial dependence of the ^-production 
operator, the FSI effect is slightly model-dependent. We estimate it by taking the ratio of 



F SD (q) = 0.968 + - ; ' 2 2 , (4) 



the Paris ^So wave function at its maximum at r = 1 fm to the corresponding plane 
wave function; numerical calculations with the Paris wave functions show that Coulomb 
effects are negligible under the conditions of the present experiment. There is a common 
enhancement factor Fss{q) for the {Ss) 2 and Ss-Sd interference terms and an analogous 
Fsd{q) for the Ss-Ds interference. They may be parametrised in a similar manner to that 
given in ref. j!2j 

151.7 

F ss (q) = 0-440+ 

1 + q 2 /a z 

11.5 
1 + q 2 /a 

with a = -0.053 fm -1 . 

We make the assumption that, apart from the FSI, the amplitudes A$ s , A$d and Ad s 
appearing in Eq. (j^J are constant. The resulting form of the differential cross section is 

N ( k 2 q 2 \ 

da = — F ss (q) + a F ss (q) cos 2 9 V + b F SD (q) cos 2 6 PP dLips , (5) 

p \ fjLrrip fim p J 

where iV is a normalization constant and dLips is the invariant three-body phase space dis- 
tribution. Typical momentum factors \xm v are explicitly shown to leave two dimensionless 
constants a and b to be determined from the shapes of the angular and energy distribu- 
tions. It must be stressed that the functional form of Eq. (jSJ does not depend upon the 
details of a specific dynamical model but on the assumption of constant amplitudes. 

To extract differential distributions from the data, we fitted the Monte Carlo simula- 
tions, with events weighted according to Eq. (4), to the experimental results. Following 
this procedure at Q = 37 MeV, a combined fit to the 6 n and Q pp angular distributions yields 
the parameter values a = — 15 ± 5 and b = 22 ± 7, which are correlated, and the overall 
normalisation constant. In addition to the given statistical errors, there are systematic 
errors, arising mainly from uncertainties in the geometrical alignment of the apparatus, 
which are of the order 10% for the a parameter and 20% for the b parameter. These values 
were then used in the Monte Carlo simulation to make the acceptance corrections needed to 
unfold the data; the resulting corrected experimental data and fitted angular distributions 
are shown in Figs, la and 2a. It should be noted that the a and b parameters were not 
deduced by fitting to the corrected data shown in the figures. 



The distribution in the rj kinetic energy T v , shown in Fig. 3a, is shifted towards higher 
energies with respect to phase space. This is a direct consequence of the strongly attractive 
pp FSI, which enhances events where the T] recoils against the two protons which have low 
excitation energy. The dashed curve, which corresponds to phase space modified by the 
FSI, produces too big an effect. The solid curve is the prediction with the values of a and 
b as determined by the angular distributions. The non-vanishing of the averages of cos 2 9 
leaves k 2 and q 2 terms in Eq. (|SJ) which yield a much better fit. 

At Q — 16 MeV, S-waves are even more dominant and the angular distributions are 
broadly compatible with phase space modified by the pp FSI. Nevertheless the shapes are 
in fact marginally better reproduced with the parameter values derived from the higher 
energy data and the corresponding predictions are shown in Figs, lb and 2b. The T v data 
shown in Fig. 3b, while again demonstrating the effect of the pp FSI, are well described by 
the parametrisation. The experimental points do depend, to some extent, on the shape of 
the assumed differential distribution, but this is well within the statistical uncertainty of 
the data. The poor quality of the data for cos6 pp > 0.8 and at the upper end of the T n 
distribution reflects mainly the loss of events due to the beam pipe so that discrepancies 
with the parametrisation should not be taken seriously in these regions. 

The numerical values of the differential cross sections are given in Tables 1-3. The data 
were normalised to the values of the total pp — > pprj cross sections given in Ref. j^j, namely 
a(T p = 1293 MeV) = (2.11 ± 0.32) fib and a(T p = 1352 MeV) = (4.92 ± 0.74) fib. 

All existing theoretical models describing reproduction in proton-proton scattering are 
broadly similar ^01 , consisting of a meson exchange exciting a nucleon isobar, dom- 
inantly the Sn N*(1535), which decays into an s-wave 77-proton pair. The models differ 
mainly in their choice of mesons exchanged and, in particular, the relative importance of 7r 
and p exchange PUJ EI] . If, for simplicity, only the pion exchange term is retained then it 
follows from the expansion of the corresponding propagator and vertex function that a = 1 
and 6 = 4 ^3]. Our experimental values are significantly larger than these. 

Higher partial wave N* resonances are potentially very important for the rj angular 
distribution. Away from threshold there is evidence for significant c?-wave production of 
the 77 meson in the 7r~p — > rjn reaction jT^j and this is confirmed in the preliminary high 

5 



statistics data from the Crystal Ball collaboration ^B]- Whilst having no appreciable 
effect on the pp angular distribution, the inclusion of such <i-wave production in a one 
pion exchange model would contribute about a ~ +7 in the r\ angular distribution which, 
though of the right order of magnitude, is of opposite sign to what we have deduced from 
our data. Thus, in contrast to our findings presented in Fig. lb, such a term would favour 
production towards cos 9 rj = ±1. 

However the sign of the cos 2 9 n term is negative in the photoproduction 7p — > rjp a 
little above threshold jTTj and, using vector dominance ideas, this is likely to be true for 
pp — > rjp as well. In a pure p-exchange model, the elementary distribution would contribute 
a m —2.5 which, though too small, is of the same sign as the one apparent in our data. 
In the original estimation of this process PHI it was claimed that p-exchange should be 
stronger than 7r-exchange and that the interference between them was mainly destructive. 
If this were indeed the case, the a coefficient could be enhanced significantly because of the 
negative sign between the p and n amplitudes. Our data would support such a conclusion. 
This would also lead to the prediction that the pn —>■ drj should show a much flatter 
distribution since the p and 7r amplitudes add in this case [TU] . 

In conclusion, we have presented the first experimental evidence for non-isotropy in 
the pp —>■ ppr] angular distributions close to threshold. The signals are generally small in 
the data as compared for example to the proton-proton final state interaction, which has 
overwhelming importance. Nevertheless, a clear sign of an rj angular anisotropy has been 
found which could be the first direct indication of p-dominance in rj production. 

An economic description of all our distributions, sufficient for acceptance corrections, 
has been given in terms of the two free parameters of Eq. © by taking the amplitudes 
A Ss , A Sd and A Ds to be constant. If, for example, we assume instead that A Ss of Eq. (0) is 
constant then the cos 2 #'s in Eq. (jSJ) are replaced by Legendre polynomials P2(cos8). The 
fits to the angular distributions in Figs. 1 and 2 are very similar but, due to the vanishing 
of the angular average of P2(cos6), the results in Fig. 3 are identical to the broken curve, 
which represents just phase space and FSI. To restore the previous good agreement requires 
an additional free q 2 or k 2 term in the fitting and, when this is introduced, the results are 
essentially identical to the solid curve of Fig. 3. The shapes of the angular distributions 
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are unaffected by such a modified procedure. 

A significant improvement of the statistics would be welcome to tie down the model 
parameters in this area. Such an improvement is anticipated through the use of the WASA 
detector |18j . which is designed for the study of rare decays of the r\ meson and which 
will have an almost Att coverage of the photons from the r\ decay. The effects of higher 
partial waves should increase strongly with beam energy but to exploit this would, in our 
case, require an energy upgrade of CELSIUS or the study of quasi-free production on the 
deuteron. 

We are grateful to the TSL/ISV personnel for their continued help during the course of 
this work. Discussions with M. Gargon on the analysis of this experiment were very help- 
ful. Financial support for this experiment and its analysis was provided by the Swedish 
Natural Science Research Council, the Swedish Royal Academy of Science, the Swedish In- 
stitute, the Bundesministerium fur Bildung und Forschung (06TU886), Deutsche Forschung 
Gesellschaft (Mu 705/3 Graduiertenkolleg) , the Polish State Committee for Scientific Re- 
search, the Russian Academy of Science, and the European Science Exchange Programme. 
The data presented here are based upon the work of J. Dyring, in partial fulfillment of the 
Ph.D. requirements [B]. 
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Table 1: Differential cross section with respect to the r] cm. angle for the pp — > ppr] reaction 
at Q — 16 and 37 MeV (T p = 1296 and 1350 MeV). In addition to the statistical error, 
there is an overall systematic uncertainty of about 20%. 



Q = 16 MeV 


Q = 37 MeV 


cos 9 V 


do j dVt^ (fib j sr) 


cos 6 V 


du / dVt ri (/ib/sr) 


-0.9 


0.170 ±0.029 


-0.9 


0.234 ±0.027 


-0.7 


0.168 ±0.033 


-0.7 


0.317 ±0.040 


-0.5 


0.123 ±0.029 


-0.5 


0.392 ± 0.049 


-0.3 


0.219 ±0.039 


-0.3 


0.552 ±0.061 


-0.1 


0.186 ±0.036 


-0.1 


0.479 ± 0.057 


0.1 


0.262 ±0.044 


0.1 


0.490 ±0.057 


0.3 


0.164 ±0.032 


0.3 


0.434 ±0.060 


0.5 


0.157 ±0.031 


0.5 


0.554 ±0.069 


0.7 


0.178 ±0.034 


0.7 


0.425 ±0.061 


0.9 


0.100 ±0.024 


0.9 


0.188 ±0.038 
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Table 2: Differential cross section with respect to the proton-proton angle for the pp — > pprj 
reaction at Q — 16 and 37 MeV. In addition to the statistical error, there is an overall 
systematic uncertainty of about 20%. 



Q = 16 MeV 


Q = 37 MeV 


cos 6 pp 


da/dQ pp (/ib/sr) 


cos 6 pp 


da/dQpp (/ib/sr) 


0.05 


0.196 ±0.039 


0.05 


0.331 ±0.043 


0.15 


0.150 ±0.029 


0.15 


0.308 ±0.041 


0.25 


0.166 ±0.028 


0.25 


0.343 ± 0.043 


0.35 


0.182 ±0.026 


0.35 


0.405 ± 0.047 


0.45 


0.169 ±0.024 


0.45 


0.344 ± 0.042 


0.55 


0.190 ±0.029 


0.55 


0.355 ± 0.042 


0.65 


0.125 ±0.029 


0.65 


0.473 ± 0.050 


0.75 


0.163 ±0.043 


0.75 


0.391 ±0.050 


0.85 


0.169 ±0.084 


0.85 


0.508 ±0.070 






0.95 


0.88 ±0.19 
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Table 3: Differential cross section with respect to the 77 cm. kinetic energy for the pp — > ppr] 
reaction at Q — 16 and 37 MeV. In addition to the statistical error, there is a typical 
systematic uncertainty of about 20%, which increases at the upper end of the spectrum. 



Q = 16 MeV 


Q = 37 MeV 


T 


da/dT v (fib/MeV) 


T 


da/dT v (/xb/MeV) 


0.5 


0.046 ±0.010 


1.0 


0.055 ±0.010 


1.5 


0.074 ±0.013 


3.0 


0.110 ±0.014 


2.5 


0.089 ±0.015 


5.0 


0.115 ±0.015 


3.5 


0.114 ±0.018 


7.0 


0.142 ±0.017 


4.5 


0.114 ±0.019 


9.0 


0.164 ±0.020 


5.5 


0.146 ±0.025 


11.0 


0.175 ±0.022 


6.5 


0.173 ±0.033 


13.0 


0.160 ±0.021 


7.5 


0.191 ±0.045 


15.0 


0.244 ±0.029 


8.5 


0.059 ±0.029 


17.0 


0.183 ±0.028 


9.5 


0.237 ±0.090 


19.0 


0.155 ±0.030 


10.5 


0.196 ±0.098 


21.0 


0.216 ±0.041 






23.0 


0.214 ±0.047 






25.0 


0.212 ±0.059 






27.0 


0.160 ±0.080 



12 



Figure Captions 



Fig. 1. Differential cross section in the ^-production angle for the pp — > ppr] reaction at (a) 
Q = 37 MeV (T p = 1350 MeV) and (b) Q = 16 MeV (T p = 1295 MeV). The dashed curves 
represent the Monte Carlo predictions of phase space modified by the proton-proton final 
state interaction, whereas the solid curve includes also the angular dependence of Eq. (5) 
with a = —15 and b = 22. 

Fig. 2. Differential cross section in the proton-proton production angle for the pp — > pprj 
reaction at (a) Q = 37 MeV and (b) 16 MeV with curves as described in Fig. 1. 

Fig. 3. Distribution in rj kinetic energies from the pp — > ppr] reaction at (a) Q = 37 MeV 
and (b) 16 MeV. The short-dashed curves represent the Monte Carlo predictions of phase 
space and the long-dashed shows the influence of the proton-proton final state interaction. 
The solid curve includes also the modifications induced by the angular and momentum 
dependence of Eq. (5) with a = — 15 and b = 22. Since the experimental acceptance 
changes somewhat according to the Monte Carlo generator used, the experimental points 
would be slightly lowered if we had extracted them using a = b = 0. 
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